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A numerical method based on the transfer matrix technique is developed to calculate the lumi- 
nescence spectra of complex layered structures with photon recycling. Using this method we show 
a strong dependence of the emission spectra on the optical eigenmode structure of the device. The 
enhancement of the photon recycling and the LED external efficiency can be achieved by placing 
the active regions inside single or coupled microcavities. 



I. INTRODUCTION 



Recently, great interest has been shown in the use of light emitting diodes (LEDs) as a light source for illumination 
jjj . LEDs offer many potential advantages compared to conventional light sources due to their relatively low energy 
consumption, long lifetime and high shock resistance. 

Currently white- light LEDs use photo-excitation of phosphors to convert the blue light from an InGaN/GaN LED 
into white light. However, phosphor has a broad emission spectrum, thus, white LEDs based on this principle do not 
have the maximum possible luminous efficacy. 

The aim of this work is to investigate the feasibility and limitations of creating a white LED by integration within 
the same structure of several semiconductor layers emitting three basic colours. The recent progress in the growth 
of InGaN-based double heterostructures and quantum wells makes this alloy an ideal material for the LED active 
regions, due to the wide variety of the energy gaps in the InGaN system, covering frequencies from red to ultraviolet. 
The working regime of such a device can be achieved by electrical pumping of active layers with the widest bandgap 
(blue regions), and by making use of re-emission of the light, absorbed by all the active regions (so called photon 
recycling) . 

II. MODELLING OF PHOTON RECYCLING 



We start by writing an expression for the intensity of spontaneous emission from a quantum well (QW) into a bulk 
dielectric material. The number of photons with energy within the interval [huj, Hlu + d(huSj\, which are emitted 
from the surface dS into the solid angle dflo during the time interval dt is given by 



dN = W dt d(hu) dS dil 
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Here E 2D is the energy gap between the electron and hole quantized levels in the QW, gin is the reduced two- 
dimensional density of states, f e and fh are the electron and hole occupation probabilities, e is the dielectric constant 
of the material containing the well, mo is the free electron mass, and Y accounts for the interband relaxation and other 
broadening mechanisms. The squared momentum matrix element p 2 v is given in the effective mass approximation for 
deep QWs by @ 



Pi 



for TE modes, and by 
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for TM modes. Here 9 is the angle between the plane- wave propagation direction and Z-axis (normal to the QW 
plane), and 7 = (E2D ~~ E g )/(huj — Eg), where E g is the bandgap energy of the QW material, and po is the interband 
momentum matrix element. 
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We will use the transfer matrix technique in the plane waves basis || to analyse the optical properties of layered 
structures. The plain-wave mode interaction with the QW is described by a QW transfer matrix. If the quantum 
well width is much smaller than the light wavelength, the transfer matrix has the form: 

M=( 1 + Y Y 
\ -Y 1 - Y 

Here Y is defined by the two-dimensional QW optical susceptibility X2D- 

Y = i-^ X 2D, (2) 
where k is the wavevector of light, fco — lu/c, and X2D for a single-subband QW is given by j| 

= K^ P l v92D l-h-h d£ 

X2D ml J (E 2D +e) 2 (E 2D +e-hLu-iT) 

To calculate the rate of photon extraction from a complex structure we consider the interference of the all possible 
processes, resulting in light emission out of the structure. These calculations require a knowledge of the amplitude 
transmission and reflection coefficients t; and 77 for the structure part on the left of the QW, and similar coefficients t r 
and r r for the structure part on the right of the QW. We also need to know the transmission and reflection coefficients 
tQw an d fQw f° r the quantum well itself. Each of these coefficients can be obtained from the corresponding transfer 
matrix. 

Let us derive, for example, the power, emitted from the right side of the structure. A photon emitted inside the 
structure to the right can be transmitted directly to the outside medium, or it can be consecutively reflected from 
the right and left sides of the structure and finally will be transmitted outside and so on. The outgoing electric field, 
resulting from all these processes is given by the sum: 

— - — t/p I 7 f T" ^ t<r | T^yT*^ ^"f^*/ I •■ — ~ - 



where the star in r* indicates, that this coefficient includes the reflection from the emitting QW. 

Similarly, photons emitted inside the structure to the left can undergo multiple reflections and eventually escape 
from the structure to the right. These processes give a second part of the external field: 

c ritQ W t r 



(1 - r*r r )(l - r QW ri) 
Thus, we obtain the expression for the emission from the right side of the structure: 
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where the ratio y^k zr /y/eok z o accounts for the change in solid angle that is due to refraction for plane waves. The 
indices r and are related to the outside medium and the layer containing the QW, respectively. To obtain the total 
density of the external light intensity we have to sum over all QWs and integrate Eq. (|J) over the external solid angle 

n e . 

If we neglect the reflection from the QW by substituting tQw = 1 and row = into Eq.(||), the formula for 
extraction becomes analogous to one obtained using the source-term method |3| . 

For the quantitative description of the recycling process, we have to calculate the rate of absorption by the quantum 
well QW a of the spontaneous emission from the other well, QW C . Thus, we need to know the induced electric field at 
the position of QW a . The power flux balance shows, that the rate of absorption of emitted photons by the unit area 
of surface of QW a is given by: 

W a - ~V~ek z \E\ 2 Re(Y QWa )W 0QWc , (4) 

where E is the complex amplitude of the field at the QW a location, Yqw is defined by Eq.(||), and Woqw s is defined 
by Eq.(||). The induced field can be calculated in a similar fashion as it was done for extraction: 



(1 - r*r* r )(l - r* al r ar ) \ 1 - r er r QWt 



2 



where r a i (r ar ) is a reflection coefficient for the part of structure to the left (right) of the well QW a , r er is a reflection 
coefficient for all the layers to the right of the well QW C , and the coefficients r and t correspond to the part of the 
structure between QW C and QW a . Here we assume that QW a is on the left of QW . The formula for the opposite 
case is similar. 

One of the channels for the photon to escape from the recycling process is to be absorbed in a metallic mirror. 
Placing such a mirror onto the left side of the structure, and denoting the reflection coefficient of the left part of the 
structure, excluding the mirror, as n, the reflection coefficient from the mirror as r m , the reflection coefficient for the 
wave incident from the mirror as r s , and the transmission coefficient from the QW to the mirror as t, we obtain the 
following expression for the number of absorbed photons: 
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We restrict our consideration by relatively low pumping levels, which are typical for the diode operation regime. 
Thus, the expression for the electron density jj] in the one-subband QW can be simplified (temperature is measured 
in energy units): 
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and the occupation probabilities can be expressed as: 
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where m e (rrihh) and /i e (uh) are the electron (heavy hole) effective mass and quasi- Fermi-level, respectively. Under 
the low-pumping assumption we can rewrite Wo in the form: 
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The total rate of light emission into the external medium can be obtained by substituting Eq.(^) in Eq.(|l|) and 
integrating Eq.(|l]) over the external solid angle and photon energies. As a result the rate of radiative recombination 
depends on the carrier densities as R ex t — E np. 

Due to the large photon energy separation between different colours, the absorption rate for the short- wavelength 
light is independent of the carrier density in the narrow-gap QW. For the light emitted from QW C the rate of absorption 
in QW a is given by the expression R a bs = An e p e , where A is calculated by integrating Eq.(Q) over the total solid 
angle of emission and over the photon energies. Here we assume that all QWs are embedded in layers of dielectric 
material with the refractive index equal to the highest one in the real structure. This trick allows as to handle the 
interaction of active layer with the evanescent wave. However, introducing several sufficiently thin layers with high 
dielectric constant does not alter the optical properties of the structure. 

A similar relation holds for the rate of absorption in a metallic mirror: R m = M np. Coefficient M is obtained by 
substituting Eq.([|) in Eq.lJHj) and integrating the result over the solid angle and photon energies. Note, that due to 
charge-neutrality the electron and hole densities in each QW are equal to each other, n = p. 

The steady-state carrier densities are given by the balance between generation of electron-hole pairs in the QWs 
and their recombination, both radiative and non-radiative. The generation processes include electric current pumping 
of the blue QWs and the re-absorption of emitted light throughout the structure. The recombination output goes to 
external emission and internal losses, which we treat as absorption in the other QWs and non-radiative recombination 
in the given QW. In our approach this leads to the following equation for each QW: 

Ri = Ein 2 + A t n 2 + Minf + m/n, 

where Ri is the pumping rate and Tj is the non-radiative recombination time. Assuming that Tj does not depend on 
the carrier density, we get: 
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Then the external efheiency for each QW can be obtained as 

= Erf 
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III. RESULTS AND DISCUSSION 



We used the numerical method described above to investigate a number of different types of structures, ft was 
revealed, that the emission spectra and external efficiencies depend drastically on the active layer positions and on the 
mode structure of the device. Figure 1 shows the calculated spontaneous emission for a structure, which represents 
GaN microcavity containing three QWs and covered on one side by a metallic mirror. The width of the GaN layer and 
the energies of interband transitions in red, green and blue QWs are chosen in such a way that the structure can be 
regarded as the 5A/2 resonator for the normally propagating red-light waves, 6A/2 resonator for the green-light waves 
and 7A/2 resonator for the blue-light waves. We placed each QW in an antinode of a resonant mode, associated with 
the QW colour. Our calculations were performed for different times of non-radiative recombination ranging from 1 ns 
to 10 ns, and we assumed the carrier density in blue QW to remain constant. One can see, that if the internal losses 
are high, the recycling efficiency is low and a blue light only is emitted. When the internal efficiency increases, the 
emission from the optically pumped QWs becomes comparable with the blue-light intensity. The green-light intensity 
usually remains smaller than both blue and red because of the strong re-absorption in the red active region. However, 
this is beneficial for the white-light generation, because of the high sensitivity of the human eye in the green region 
of the spectrum. If high non-radiative loses are present, the recycling process can be enhanced by introducing more 
red and green QWs in the structure and by building a Bragg reflector for the blue wavelength. A possible way to 
enhance the external efficiencies of all three colours is to place the active regions into coupled microcavities. 
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FIG. 1. Spontaneous emission spectra from a microcavity, containing three quantum wells. The non-radiative recombination 
times are 1 ns (solid line), 5 ns (dashed line) and 10 ns (dot-dashed line). The carrier density in the blue QW is the same for 
all three curves. 
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